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Observationally, both the 3.4 /im aliphatic hydrocarbon C-H stretching absorption feature and the 9.7 /im amor- 
phous silicate Si-O stretching absorption feature show considerable variations from the local diffuse interstellar 
medium (ISM) to Galactic center (GC): both the ratio of the visual extinction (Ay) to the 9.7 /im Si-O optical 
depth (Arg.7 Mm ) and the ratio of Ay to the 3.4 C-H optical depth (Ar3.4^m) of the solar neighborhood local 
diffuse ISM are about twice as much as that of the GC. In this work, we try to explain these variations in terms of a 
porous dust model consisting of a mixture of amorphous silicate, carbonaceous organic refractory dust (as well as 
water ice for the GC dust). 
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1. Introduction 

The interstellar extinction law is one of the primary 
sources of information about the interstellar grain popula- 
tion, and one often obtains direct information on the compo- 
sition of interstellar dust from spectral features in extinction 
(Draine 2003). These spectral features also provide strong 
constraints on interstellar grain models. With the advent of 
ground-based and space borne infrared (IR) telescope facil- 
ities, the IR extinction continuum and absorption features 
have been receiving increasing attention and play an essen- 
tial role in recovering the intrinsic energy distribution of ce- 
lestial objects and inferring the characteristics of interstellar 
dust. 

In the interstellar extinction curve, the 2175 A bump is 
outstanding in the ultraviolet (UV), while in the IR, there are 
a number of prominent absorption features as well: (1) the 
ubiquitous 9.7 /im and 18 /im features respectively due to the 
Si-O stretching and O-Si-O bending modes of amorphous 
silicates; (2) the 3.4 /im feature due to the C-H stretching 
mode of aliphatic hydrocarbon dust, as ubiquitously present 
in the ISM of the Milky Way and external galaxies as the 
9.7 /tm and 18 /im silicate bands, except this feature is not 
seen in dense molecular clouds (see Pendleton 2004 for a 
review); (3) the 3.3 /im and 6.2 /im weak features seen in 
both local sources and GC sources (Schutte et al. 1998; Chiar 
et al. 2000), respectively due to the C-H stretching and C-C 
stretching modes of polycyclic aromatic hydrocarbon (PAH) 
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molecules; and (4) in dense clouds the 3.1 /im feature due to 
the O-H stretching mode of water ice as well as a number of 
weaker features at 4.68 /im (CO), 7.68 /im (CH 4 ), 4.28 /im, 
15.2 /im (C0 2 ), 3.54 /im, 9.75 /mi (CH 3 OH). 

The 9.7 /im silicate extinction profile varies among differ- 
ent sightlines; in particular, its optical depth Atqj^ (rel- 
ative to the visual extinction Ay) shows considerable vari- 
ations from the local diffuse ISM (LDISM) to the Galactic 
center (GC): Ay / At 9 . 7 ^ « 18.2 for LDISM differs from 
that of the GC (Ay / At 9 , 7 ^ » 8.4) by a factor of ^2.2 
(see Table 1; also see Draine 2003)Q Roche & Aitken (1985) 
argued that Ay / Arg^m varies because there are fewer car- 
bon stars in the central regions of the Galaxy and the produc- 
tion of carbon-rich dust may be substantially reduced com- 
pared with the outer Galactic disk. However, as shown in Ta- 
ble 1, the 3.4 ftm C-H feature of aliphatic hydrocarbon dust 
also exhibits a similar behavior: Ay /At^^^ k 274 for 
LDISM is higher than that of the GC (Ay / At 3 4 Mm w 146) 
by a factor of ~ 1.9 (see Table 1). If the argument of Roche 
& Aitken (1985) was valid, one would expect a much smaller 
Ay I Ar 3 . 4 Mm ratio in the LDISM than that of the GC. Sand- 
ford et al. (1995) tried to quantitatively explain this phenom- 
ena by assuming that the abundance of the C-H carrier (rel- 
ative to other dust components) gradually increases from the 



1 Cohen et al. (1989) argued that the observed 9.7 /im dip on which 
AT9.7„ m was measured may be partly contributed by PAHs (i.e. the red 
tails of the PAH 7.7 pm and 8.6 /im bands and the blue tail of the 1 1.3 pro 
band could form an "artificial" 10/im dip). But this would result in a 
smaller Arg.7 Mm for the GC and a larger Arg.7 Mm for the LDISM (since 
the PAH emission is more likely to present in the LDISM while toward the 
GC PAHs are seen in absorption), quite on the opposite. 
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Table 1. 


Observational values of A 


WAT3.4 Mm and A v 


/Arg.7 




Ay /AT3.4 Mm 


Av/Atq.7 Mm 


References 




240±40 




Sandford et al. 1991 




250±40 




Pendleton et al. 1994 


Local diffuse ISM 


333 


16.7 

18.5 
18.5±2 
19.2±0.6 


Adamson et al. 1990 
Roche & Aitken 1984 
Draine 2003 
Bowey et al. 2004 


mean for local ISM 


274 


18.2 






150±20 




Pendleton et al. 1994 


Galactic Center 


143 


8.31 

8±3 
9 


McFadzean et al. 1989 

Becklin et al. 1978 
Roche & Aitken 1985 


mean for GC 


146 


8.4 





local ISM toward the GC0 However, this requires that amor- 
phous silicate dust and aliphatic hydrocarbon dust should not 
be solely responsible for the visual extinction. If one has to 
invoke an additional dust component (most likely a popula- 
tion of carbon dust which does not show the characteristic 
3.4 fj,m feature, say, graphite) making an appreciable contri- 
bution to Ay, one would encounter a severe carbon budget 
problem (see Snow & Witt 1996). 

Along the lines of sight toward the GC, there are dense 
molecular cloudsH In cold, dense molecular clouds, inter- 
stellar dust is expected to grow through coagulation (as well 
as accreting an ice mantle) and the dust is likely to be porous 
(Jura 1980). In this work, we demonstrate that the observed 
variations of Ay / Atqj Mm and Ay / AT3.4 Mm from the LD- 
ISM to the GC could be explained in terms of composite 
porous dust. 

2. Model 

We consider a composite porous dust model consisting of 
amorphous silicate, carbon, and vacuum (in dense clouds sil- 
icate dust and carbon dust are coated with water ice). We 
take the optical constants of Draine & Lee (1984) for amor- 
phous silicate, of Li & Greenberg (1997) for carbonaceous 
organic refractory (to represent the carbon dust component), 
of Li & Greenberg (1998) for water ice. The mass densi- 



compact dust, a=0A fim 



2 They also pointed out that the C-H and Si-O carriers may be coupled, 
perhaps in the form of silicate-core organic-mantle grains. This idea is 
challenged by the nondetection of the 3.4 /an feature polarization along 
sightlines where the 9.7 fim feature polarization is detected (Adamson et al. 
1999, Chiar et al. 2006; also see Li & Greenberg 2002). It is possible that 
the 3.4 fim feature may be not produced by a carrier residing in a mantle on 
a silicate core but by very small (unaligned) grains (Chiar et al. 2006). 

3 The sightline toward the Galactic center source Sgr A* suffers about 
~ 30 mag of visual extinction (e.g. see McFadzean et al. 1989), to which 
molecular clouds may contribute as much as ~ 10 mag (Whittet et al. 1997). 



porous dust, a=0.0l fim 
porous dust, a-0. 1 fim 
porous dust, a=l .0 fim 
P = 0.8 




Fig. 1. Extinction cross sections C ex t(A) of different types of dust. The 
3.1 fim water ice O-H feature shows up in the extinction profiles of 
porous dust. 



ties of silicate dust, organic refractory dust and ice are taken 
to be p s a « 3.5 gem -3 , p carb « 1.8 gem -3 and p icc w 
1.2 gem -3 , respectively. We take the mass ratio of organic 
refractory dust to silicate dust to be m ca rb/w 8 a = 0.7 and 
the mass ratio of water ice to organic refractory dust and sili- 
cate dust to be mice/ (»ri car b + m s u) — 0.8, as inferred from 
the cosmic abundance constraints (see Appendix A of Li & 
Lunine 2003a). 

For the dust in the local diffuse ISM, we assume the dust to 
be a solid compact mixture of amorphous silicate and organic 
refractory materials with m car b/m S ii = 0.7. We take the 
dust size to be a = 0.1 /im, the typical grain size for the 
dust in the diffuse ISM (see Draine 1995). For the dust in 
the dense molecular clouds along the lines of sight toward 
the GC, we assume that silicate dust and organic refractory 
dust are equally coated with an ice layer and then form a 
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A (/im) A (/^m) 

Fig. 2. Schematic illustration of obtaining AC cx t(9.7/im), the excess 9.7 /an extinction cross section above the continuum. It is obtained by (1) 
fitting the 9.7 fj,m model profile with a Drude function and subtracting the continuum which is fitted with a six-order polynomial, (2) integrating the 
continuum-subtracted extinction profile (which is approximated by a Drude function) over wavelength, and finally (3) dividing the integrated value with 
the width of the interstellar 9.7 /im silicate absorption feature. 
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Fig. 3. Same as Figure 2 but for the 3.4 (im feature. In the right panel, the 3. 1 fim water ice O-H feature shows up in the extinction profile of porous dust. 



porous aggregate (see Li & Lunine 2003b). For porous dust, 
a key parameter is the porosity P (or fluffiness; the fractional 
volume of vacuum in a grain). We will consider a range 
of porosities. We assume all grains are spherical in shape; 
the porous grain size a is defined as the radius of the sphere 
encompassing the entire porous aggregate. In order to find 
suitable porosity P and dust size a for the dust in the dense 
molecular clouds to reproduce the observed Ay / ATg.7 pm 
and Av/^T3.4fj,m ratios toward the GC, we leave both P 
and a adjustable. 

We use Mie theory in combination with the Maxwell- 
Garnett and Bruggeman effective medium theories (Bohren 
& Huffman 1983; see eqs.7-9 of Li & Lunine [2003b] and 
Kimura et al. 2008b) to calculate the optical properties of 
composite porous grains. This approach is valid for com- 
puting the integral scattering characteristics (e.g. extinction, 



scattering, absorption cross sections, albedo and asymmetry 
parameter; see Hage & Greenberg 1990, Wolff et al. 1994). 

For illustration, we plot in Figure 1 the 1-100 (im extinc- 
tion cross sections C ox t (A) of compact dust (for the local dif- 
fuse ISM) and porous dust (for dense clouds toward the GC). 
For a given dust size, both Ay and AT3.4 Mm , AT9.7 de- 
crease with the porosity P, although the degree of decrease is 
somewhat less significant for AT3.4 Mm and AT9.7 Mm than for 
Ay (thus Ay/AT3.4 pm and Ay I A79.7 fj, m decrease moder- 
ately with the increase of P; see Table 2). This is because 
with the same size, a porous grain contains less dust material 
than its solid counterpart. The effective dielectric functions 
of porous dust are reduced. In the IR the extinction for dust 
of ~ 0.1 /Ltm is dominated by absorption (i.e. in the Rayleigh 
regime) and is roughly proportional to the imaginary parts of 
the dielectric functions (see Li 2008). Therefore porous dust 
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produces smaller Ar3.4 Mm and Arg.7 Mm than its solid coun- 
terpart of the same size. In the optical, both absorption and 
scattering are important. The introduction of vaccum leads 
to a reduction of the dielectric functions of the dust which 
will decrease both the absorption and scattering. 

For models consisting of single-sized dust, the ratio of 
Ay to the optical depth of the 9.7 /im silicate feature is sim- 
ply A v /Ar 9 . 7 ^ m « 1.086 C ext (V)/AC ext (9.7 fim), where 
C ex t(V) is the extinction cross section at V band (A = 
5500 A), AC cx t(9.7 /im) is the excess extinction cross sec- 
tion of the 9.7 /im feature above the continuum, and the fac- 
tor "1.086" arises from the conversion of extinction (in mag- 
nitude) to optical depth. We obtain AC 0X t(9.7/im) by in- 
tegrating the 9.7 /im model extinction profile (with the con- 
tinuum subtracted) over wavelength (see Fig. 2 for illustra- 
tion) and then dividing the integrated value with the width 
of the interstellar 9.7 /im silicate absorption feature. The 
same procedure is applied to the 3.4 /im feature to calculate 
AC 0X t (3.4 /im) so as to obtain Ay / AT3.4 Mm (see Fig. 3 for 
illustration). 

3. Results and Discussion 

We calculate Ay / AT3.4 and Ay / AT9.7 for vari- 
ous dust models with a range of porosities and dust sizes. 
The results are tabulated in Table 2. For a given dust size 
a = 0.1 /xm, porous dust results in smaller Ay / AT3.4 fj, m 
and Ay /AT9.7 jum values than compact dust (see §2). This 
is encouraging that porous dust which is likely present in 
the dense molecular clouds along the sightlines toward the 
GC indeed is on the right track to account for the ob- 
served Ay j At 3.4 and Ay / Arg.7 Mm variations from the 
local diffuse ISM to the GC. More specifically, from Ta- 
ble 2 we see that the dust with P ~ 0.8 — 0.9 and a ~ 
0.5 — 1 /im can approximately explain the observed varia- 
tions of Ay /Ar3.4 pm and Ay /Arg.7 Mm (by a factor of ~ 2) 
from the local diffuse ISM to the GC. Both high porosities 
(P ~ 0.8 — 09) and large sizes (a > 0.5 /im) are required 
for the GC dust to account for the lower Ay / At$ .4^ and 
Ay/Ar 9 .7 Mm ratios. 

In Figure 4, we show Ay / A73.4 Mm and Ay /Arg.7 Mm as a 
function of dust size for P = 0.8, 0.9, 0.95. It is clearly seen 
that for a given porosity, the variation of Ay / A73.4 Mm with 
dust size exhibits a tendency similar to that of Ay /Arg.7 Mm . 
This suggests that with a dust size distribution taken into ac- 
count, we would still maintain the variation tendency. For 
small, highly porous grains (a < 0.05 /im), Ay / A 73. 4^ 
and Ay / Arg.7 Mm are nearly independent of size since they 
are more or less in the Rayleigh regime in the optical-IR and 
therefore Ay /V dust , /Vdust and AT 9 .7 Mm /Vd ust 

are independent of grain size (where Vdust is the dust vol- 
ume). The visual extinction Ay reaches its maximum at 
a ~ A/ [27T a (n — 1)] (where n is the real part of the refrac- 
tive index of the porous dust at wavelength A; see Li 2008). 
At even large sizes, while Ay reaches a constant (i.e. in 
the geometric optics regime) A73.4 ^ and Arg.7 /im increase 
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Fig. 4. Variations of Ay / AT3.4 Mm and Ay / AT9.7 ^ m with dust size. 



with a (till they reach their respective peak values). This ex- 
plains why Ay /A73.4 ^ m and Ay j Arg.7 Mm decrease with a 
after they reach their peak values. 

For Ay j AT3Afj.n1, out model with m car t, /?rj s ii = 0.7 
(and P — 0.8, a = 0.5 — 1.5 /im) is consistent with 
the observed factor-of-two variations in the local ISM and 
toward the GC (see Tables 1,2). However, the model 
values of Ay / At 9 7 Mm for both the diffuse ISM dust 
(Ay / Arg.r fj, m ~ 38.2) and the GC dust (A v /At 9 . 7 ^ w 
16.3) are higher by a factor of ~ 1.5-2 than that observed 
in the local diffuse ISM (Ay /Ar 9 7Mm » 18.2) and the GC 
(Ay I Arg.7 ium ~ 8.4). This discrepancy may result from the 
underestimation of the silicate mass fraction. 

With an increased silicate mass fraction, say, 
TOcarb/TO s ii = 0.5 which is consistent with the in situ 
measurements of comet Halley (Jessberger & Kissel 1991) 
and widely adopted in cometary dust modeling (Greenberg 
1998; Greenberg & Li 1999; Kimura et al. 2006, 2008a; 
Kolokolova et al. 2004, Kolokolova & Kimura 2008; 
Mann et al. 2006), we obtain Ay / AT^A^m ~ 252 and 
Ay j Ar 9 .7 « 27.1 for the local ISM (assuming compact 
dust), and Ay / At 3A ^ m ps 154 and ^4y/Arg.7 Mm « 11.3 
for the GC (assuming porous dust). These values are 
closer to that observed. It is expected that with a smaller 
TOcarb/w s ii (i.e. a larger silicate mass fraction), one would 
obtain a smaller Ay / Arg.7 nm while Ay / At^a Mm does not 
change much. Thus the observed variations of Ay j AT3 .4 Mm 
and ^4y/Arg.7 Mm from the local ISM to the GC could 
be explained. It is worth noting that, based on a detailed 
analysis of the GC 5-8 /im absorption spectra obtained 
from the Kuiper Airborne Observatory, Tielens et al. (1996) 
argued that silicate dust may contribute as much as 60% 
of the interstellar dust volume. This would translate to 
TOcarb/rTisii « 0.34 if we assume that the remaining 40% 
of the interstellar dust volume is all from the 3.4 /im C-H 
feature carrier (which is indeed a very generous assumption). 
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Table 2. A v 


/Ar 3 .4 M m and .4 


V//Ar9.7,j m 


calculated for various dust models 




Dust 


Porosity 


size 


Av/At 3 A 




FWHM 


Av/At().7 M m 


FWHM 


Model 


P 


w 






(^m) 




(H 


diffuse ISM Dust 





0.1 


243 




0.14 


38.2 


2.18 




0.1 


0.1 


216 




0.14 


32.3 


2.19 




0.2 


0.1 


1 no 

192 




0.14 


27.3 


2.19 






0.1 


80 




0.14 


5.9 


2.18 






0.5 


150 




0.15 


13.3 


2.18 




0.8 


1.0 


186 




0.17 


16.3 


2.18 




1.5 


127 




0.19 


14.8 


2.21 






2.0 


102 




0.19 


11.5 


2.26 


GC Dust 




3.0 


55 




0.20 


6.1 


2.34 




0.1 


72 




0.14 


5.1 


2.08 






0.5 


126 




0.14 


9.3 


2.08 




0.9 


1.0 


165 




0.15 


12.8 


2.09 




1.5 


193 




0.16 


15 


2.11 






2.0 


198 




0.15 


15.6 


2.11 






3.0 


174 




0.16 


14.1 


2.16 



Diffuse ISM dust is taken to be a mixture of amorphous silicate and organic refractory substance. 



Admittedly, the proposed explanation is oversimplifed. In 
the future we will consider more realistic models in which 
more dust species (e.g. hydrogenated amorphous carbon 
with a range of C/H ratios), the distribution of dust along the 
line of sight toward the GC (e.g. see Sandford et al. 1995), 
a distribution of dust sizes, and the possible porous nature of 
the diffuse ISM dust (e.g. see Mathis & Whiffen 1989) will 
be considered. 
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